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Abstract. - We study a new statistics of wave functions in several chaotic and disordered 
systems: the random matrix model, band random matrix model, the Lipkin model, chaotic 
quantum billiard and the ID tight-binding model. Both numerical and analytical results show 
that the distribution function of a generalized Riccati variable, defined as the ratio of components 
of eigenfunctions on basis states coupled by perturbation, is universal, and has the form of 
Lorentzian distribution. 



Statistics of eigenfunctions in disordered and chaotic systems are of great interest in many 
branches of physics such as condensed matter physics, nuclear physics, chemical physics, and 
in particular, in the new developing field of quantum chaosjlj |2| [|. However, although eigen- 
functions contain more information than eigenenergies, much less has been studied compared 
with the eigenenergies. 

The mostly studied property is the statistics of wave function amplitude |^(r)| 2 an d two 
point correlation function (^(ri)^^)). The former one has been found to be the Porter- 
Thomas distribution for time-reversal invariant systems whose classical counterpart is chaotic. 
This is predicted by the Random Matrix Theory |Q ^, and confirmed numerically ^ and 
experimentally Jio| . The later one, which is proportional to Jo(kd) (Jo is the Bessel function of 
the zeroth order, k the wave vector and d the distance between the two points), was proposed by 
BerryQ for a chaotic billiard within an assumption that the wave function is a superposition of 
plane waves with random coefficients. The two point correlation function has been extended to 
large separation by Hortikar and SrednickiQ recently. A beautiful review for the correlations 
of wave functions in disordered systems has been given recently by Mirlin |1 if . 
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In this Letter, we would like study another interesting statistics, namely the statistics of 
the Riccati variable of wave functions (see, e.g., p2[). The Riccati variable has been studied 
in one-dimensional disordered system, where it is related to the reflection coefficient of waves. 
As we shall see later, the Riccati variable can be easily related to the two-point correlation 
function, but it contains more information than the correlation function. Furthermore, the 
Riccati variable gives information about local fluctuation properties of wave functions, and the 
form of its distribution is related to properties of the perturbation of the physical models. To 
seek a general universality of this quantity among different quantum systems, our study will 
be extended to a wide range of random matrix model, disordered and chaotic systems such as 
the Random Matrix Model (RMM)Jff], Band Random Matrix Model (BRMM)JIJ], the Lipkin 



model [y_5| with chaotic behavior in the classical limit, a chaotic quantum billiard model) 16 
and a ID tight-binding model in the case of weak disorder fl7f. 

In order to study the chaotic systems, the Riccati variable for the ID tight-binding model 
is first extended to a series types of Riccati variable, namely, type I, type II and so on. Then, 
the statistics for type I Riccati variable is found analytically to be of a Lorentzian form for the 
RMM. Based on this fact, we conjecture that the distribution function of the Riccati variable 
of type I is universal for all chaotic and disordered systems. Our numerical results on several 
such models support the conjecture. Furthermore, the form of the distribution of the other 
types of Riccati variable gives information on properties of chaotic eigenfunctions which can 
not be supplied by, e.g., the statistics of intensity of wave functions. 

Let's first consider a general case in which the Hamiltonians can be written in the form 

H = H + XV, (I) 

where A is a parameter for adjusting the strength of the perturbation. Eigenvectors of Ho and 
H will be denoted by \i) and \a) in what follows, respectively, i.e., 

H \i) = E° i \i), H\a)=E a \a). (2) 

Usually the perturbation V can be written in the form of V = J2t ^ f with Vt = Ut + U} and 
Ut coupling a basis state to one of the other basis states only. In this Letter, we choose the 
strength of Vt comparable with each other. The component of an eigenstate \a) on a basis 
state |i) will be denoted by C a i(= (i\ct)). For simplicity, we consider the case that both H 
and Hq are time-reversal invariant only, so that C a i are real. 

In the ID tight-binding model the so-called Riccati variable is defined as the ratio of 
components of eigenfunctions of the nearest sites, where the off-diagonal elements of the 
Hamiltonian are non-zero for the nearest sits only. A natural extension of the concept to 
the general case of H in Eq. ([[]) is to define the Riccati variable as p a ,(ij) = C a i/C a j for |i) 
and \ j) coupled by a Vt ((i|Vf|j) ^ 0). Note that the label j is determined by the label i due to 
V±- The distribution of p a ,(i,j) f°r au possible i of an eigenfunction, or of eigenfunctions with 
eigenenergies in some energy range, will be denoted by (the subscript E will be omitted 

for brevity in what follows). 

Furthermore, the ratio p a ^i.j) for \i) and \j) not coupled by V is also of interest here. 
When (i|F|i) = but (i|V 2 |j) ^ 0, the ratio p a ,(i,j) will be called type II Riccati variable. 
The Riccati variable mentioned above for \i) and \j) coupled by V, is type I. For \i) and \j) 
satisfying (i| V"]j) ^ but (ilV 1 " 1 ]^ = 0, the ratio P a ,U,j) wm be called type n Riccati 
variable. The distribution of type n Riccati variable is also denoted by f(p). Expressing 
C a iC a j as p a ,{i,j)Caji one can see that for chaotic eigenfunctions, whose components and 
ratios of components can be regarded as statistically independent, the average value of p a uj\ 
is directly related to the two-point correlation function of a fixed distance by the relation 
(Pa,(i,i)><C*i> = (C*iC aj ) where (C 2 aj ) = 1/N. 
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Random Matrix Model. - We start our discussion with the RMM|13||. Choosing A = 1 
in Eq. ([j]) and taking the matrix elements of Ho and V in the Hq representation to be real 
random numbers with Gaussian distribution {(Hf k ) — 1), one has the Gaussian orthogonal 
ensemble (GOE). The joint distribution of two arbitrary components Ci and Cj for an arbitrary 
eigenfunction is (subscript a omitted for brevity) 

p(r m - r((JV + i)/2) 2 r2){N - 3)/2 m 

(see ref. Q). Then, the distribution f(p) can be readily obtained 

f G oE(p) = J dddCAp - §:)P(Ci, = (4) 

which is a Lorentzian distribution centered at p = 0. Notice that the ratio of any two 
components of eigenfunctions of the GOE is a type I Riccati variable. 

It is then natural to conjecture that the distribution of type I Riccati variable for chaotic 
eigenfunctions also has a Lorentzian form. As will be shown below numerically for other 
models, the distribution of the Riccati variable of type I for chaotic eigenfunctions is in fact of 
a more general Lorentzian form, denoted by fh{p) in order to be distinguished from /goe(p), 
which may center at p ^ 0, 

fL{p)= (p (5) 

where a — yl — s 2 and s = (p) which is proportional to the value of two point correlation 
function for a fixed \i — j\. In fact, fh{p) can be obtained under a condition more general 
than the RMM. Consider a distribution f(p) for p and the distribution F(q) for q = 1/p. If 
F(q) has the same form as f(p), i.e., F(q) = f(q), and l/f(p) can be expanded as a Taylor 
series l//(p) = do + dip + dip 2 + • ■ ■ , then f(p) must has the form of /l(p). This can be 
proved easily by making use of the relation between F(q) and f(p) required by q = 1/p, i.e., 
F(q) = f(l/q)/q\ 

Lipkin Model. - The first dynamical model employed in the Letter is a three-orbital 
schematic shell model, called the Lipkin model, which is bounded and conservative. In this 
model there are totally f2 particles distributed in three orbitals. Both the quantum and 
the corresponding classical dynamics of the model have already been known (see, e.g., Refs. 
@ H). The Hamiltonian used in Ref. has the form of Eq. (0) with 

4 

Ho = eiK n + e 2 K 22 , V = ^/i t F t , (6) 

t=i 

where 

Vi = KioKw + KoiKoi , V 2 = K 20 K 2 o + K 02 K 02 , 

V 3 = K 21 K 20 + K Q2 K 12 , V A = K 12 K W + K Q iK 21 . (7) 

Here e; and fit are parameters given in Jl9| , for which the system is almost chaotic when A = 2. 
The operators K o,Kn and K 22 are particle number operators of the orbitals 0, 1 and 2, and 
K rs for r s are particle raising and lowering operators, respectively. The eigenstates of H a 
are \mn) — A(m, n)K"oK2o\00), where A(m,n) is the normalization coefficient, m and n are 
particle numbers of the orbitals 1 and 2, respectively. The total particle number is conserved, 
and as a result the particle number of the orbital is £1 — m — n. In our calculations we take 
SI = 40 and the dimension of the Hilbert space is JV = (fi + l)(fi + 2)/2 = 861. 
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The basis states \mn) can also be labeled in energy order as In ^-representation, the 
Hamiltonian matrix is banded [[Ts(] . For convenience, we can put the two labels together with 
rrii and rij indicating particle numbers of the orbital 1 and 2, respectively. Eq. (R) tells that 
the perturbation couples only four kinds of basis states with definite values of Am — rrii — rrij 
and An = n^ — rij, e.g., Am = ±2, An = for V\. That is, there are totally four kinds of 
Riccati variable of type one. 

When A = 2, in the middle energy region the nearest-level-spacing distribution of the 
Lipkin model is close to the Wigner distribution and the underlying classical dynamics is 
almost chaotic (except a few quite small regular islands). For this perturbation strength, it 
has been found that the distributions of all the four kinds of Riccati variable of type I are 
quite close to the distribution As an example, in Fig. 1(a) we present the distribution 

f(p) (dots, in the logarithm scale) for p — C a i/C a j with \i) and \j) coupled by V\. In order to 
have a good statistics we have diagonalized 51 Hamiltonian matrices with |A — 2| < 0.025 and 
then put data obtained for eigcnfunctions of a = 430 — 450 together. The solid line in Fig. 
1(a) is the best fitting curve of the Lorentzian form Jl{p)- The central parts of the f{p) and 
the fitting curve are presented in the inset of Fig. 1(a) in the ordinary scale. We would like 
to mention that the center of the distribution f(p) in Fig. 1(a) is not at p = 0. 

On the other hand, for Riccati variable of type n with n > 2 (also for A = 2), it has been 
found that the distribution f(j>) deviates from the Lorentzian form /l(p), and the larger n 
the larger the deviation is. In Fig. 1(b) an example is given for the distribution of type V 
Riccati variable with TOj — mj — ±7 and n, — rij = 0, where deviation of f(p) (dots) from the 
Lorentzian fitting curve (solid line) is quite obvious. However, when the perturbation strength 
A is increased further, the deviation from Jl(p) will become smaller. In fact, for each type 
of Riccati variable, so long as the value of A is large enough, the distribution function will 
become Lorentzian. For example, when A = 6 the form of f(p) for type II Riccati variable can 
be fitted as well as the f(p) in Fig. 1(a) by the Lorentzian distribution Jl{p)- 

Finally, we would like to mention that if the system is not chaotic the distribution f(p) for 
all types of Riccati variable deviates from the Lorentzian distribution The deviation 

increases with decreasing A. When A approaches to zero, namely, the system becomes near 
integrable, the distribution f(p) will approach to a S function. 

Band Random Matrix Model. - Since the difference between the Lipkin model and the RMM 
lies in the band structure of the Hamiltonian matrix of the Lipkin model in Hq representation, 
we have also studied the distribution f(p) for the Band Random Matrix Model (BRMM). 
For this model we take A — 1 in Eq. (|lj) and Hik being real random numbers with Gaussian 
distribution ((Hf k ) — 1) for \i — k\ < 6 and zero otherwise. When b > 1, despite of whether the 
eigenfunctions are localized or extended, numerically it has been found that the distribution 
f(p) for Riccati variable of type I j\ < b) is very close to /gos(p) (due to the randomness 
of the center of f(p) is at p = 0); while for the distribution of Riccati variable of type n 
with n > 1, deviation of the f(j>) from /gO-e(p) has been found which enlarges gradually as n 
increases. An example is given in Fig. 1(c) for i — j = ±46, i.e. type IV in which the difference 
between the distribution f(p) (dots) and the fitting curve of the Lorentzian form (solid line) 
is quite clear. When A is increased, similar to the case of the Lipkin model, f(p) distributions 
for Riccati variable of type n with n > 1 will also approach to the Lorentzian form of fcOE(p)- 

In order to compare with the result for the Lipkin model that centers of f(p) are not 
necessarily at p = 0, we have tried to add a constant to all the non-zero off-diagonal elements 
of the BRMM. Then, it has been found that the peak of f(p) could indeed be shifted (see 
Fig. 1(d)) to s(t^ 0). The value of s has been found numerically in a good agreement with 
two-point correlation function. 

Chaotic Quantum Billiard. - Chaotic quantum billiards have been widely studied as a 
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Fig. 1. - (a) The distribution function f(p) (dots) of type I Riccati variable for the Lipkin model in 
chaotic region, (b) Same as (a) but for type V. (c) The f(p) of type IV Riccati variable for the Band 
Random Matrix model, i.e. i — j = 20 band width 6 = 5. (d) The f(p) (dots) obtained by adding a 
constant to all the non-zero off-diagonal matrix elements of BRMM. The solid curves in (a)-(d) are 
best fitting curves of the Lorentzian form. 

prototype of quantum chaos. We take the Bunimovich stadium|l6| as an example. A chaotic 
odd-odd wave function of a stadium with radius R = 1 and a = 1 has been calculated by 
plane wave decomposition method|8). The wave function is shown in the inset of Fig. 2(a). 
For this wave function the statistics of the intensity 2 is a Porter- Thomas distribution. 
The distribution f(p) for p = ip(x,y)/ip(x,y + 5y) with a fixed 8y larger than the de Broglie 
wavelength is shown in the figure, which is found to be a very good Lorentzian. The best 
fit by Eq. ([)]) gives rise to s = (p) = 0.16. As was pointed out before, s is related 
to the two point correlation function at a fixed 5y. According to Berry's conjecture^, 
(ijj(x,y)il}(x,y + 5y)) = Jo{kSy), where A (= l + 7r/4) is the area of the billiard, and k5y — 10n 
here. Thus, we have (p) = AJo(k8y) — 0.18. This value is very close to the one obtained 
above. 

ID Tight Binding Model. - Finally, although statistical properties of spectra and distribu- 
tion of components of eigenfunctions for classically chaotic quantum systems are quite different 
from those of the ID tight binding model (Anderson model) in the case of weak disorder, it 
has been found that for the tight-binding model the distribution of Riccati variable of type 
I is also close to the Lorentzian form. For this model we take = 5i^+i + <5;,fc-i and E® 
being random numbers with flat distribution in the region [—1/2, 1/2]. An example is given 
in Fig. 2(b) for the case of A = 1.0. The agreement between the f(p) distribution (dots) and 
the fitting curve of the Lorentzian form (solid lines) is quite obvious. 

In conclusion, for all different chaotic and disordered models studied in this Letter, we have 
shown that the distribution f(p) for type I Riccati variable is universal and has the form of 
Lorentzian distribution. On the other hand, the distribution function for Riccati variable 
of type n with n > 1, although approaching to the Lorentzian form when perturbation 
is strong enough, shows different features in different stages of perturbation strength even 
though the systems are already chaotic. Therefore, statistics of Riccati variables of different 
types reflects statistical properties of eigenfunctions that can not be reflected by the statistics 
of wave function amplitude. Furthermore, the properties of Riccati variable may supply an 
alternative measure to probe information for coupling structure (the Hamiltonian structure) 
for real systems. 
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